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ABSTRACT 



The design, testing and analysis of a dual frequency 
system to detect and determine the resonance frequency (and 
hence the size) of bubbles ranging from 100 to 7 jum radius is 
reported. The resonance frequencies were compared to 
estimates based on the rise time of the bubbles. In general 
these comparisons agreed to within five percent. Although the 
system is not ideal for field measurements, it identifies 
important requirements concerning the size of the sample 
volume, the frequency and amplitude ranges of the sound fields 
and signal processing techinques to make an efective system. 
Bubbles are distinguished from non-gaseous particles by the 
nature of their nonlinear response to the dual sound field. 
The system is versatile and can be modified to suit many 
research purposes. 
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I . INTRODUCTION 

Bubbles in the ocean create a wide variety of effects on 
sound propagation including absorption and scattering of 
acoustic energy. In order to develop an accurate model of 
these effects, it is important to determine the number and 
size of bubbles in the regions of interest. This thesis 
reports on an implementation of the dual frequency method to 
detect bubbles and determine their resonance frequency. 

This implementation was conceived with the following goals 
in mind: 

* Detect bubbles with resonance frequencies from 30 kHz to 
400 kHz using the dual frequency method published by 
Newhouse and Shankar [Ref. 1] and investigated by Hampton 
[Ref. 2]. 

* Measure the resonance frequency of bubbles as they are 
detected . 

* Distinguish between bubbles and non-gaseous particles. 

* Determine the parameters required for a high probability 
of detection for bubbles in the sample volume. 

The dual frequency method is explained in the theory 
section. Supporting theories regarding bubble resonance, 
visual determination of bubble size and nonlinear mixing 
properties of solids are included. The problem approach 
section details the experimental apparatus and tests. The 
results of the measurements are then presented and discussed. 
This thesis concludes with a summary and recommendations for 
future work in this area. 
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II 



THEORY 



The dual frequency method of bubble detection and sizing 
relies on the nonlinear properties of bubbles undergoing large 
amplitude oscillations. Theories concerning the 
characteristics of bubble resonance and nonlinear signal 
generation by bubbles in a dual frequency sound field are 
presented to clarify the dual frequency method. Also, because 
one of the goals of this project is to distinguish between 
bubbles and solid particles, it is necessary to understand how 
sound is scattered by solid particles exposed to a dual 
frequency sound field. A discussion of radiation pressure, 
streaming and rectified diffusion is included since these 
mechanisms limit the allowable pressure amplitudes. A section 
on bubble rise time sizing covers the method used to verify 
the resonance frequency of certain bubbles by an alternate 
means . 

A. BUBBLE RESONANCE 

The fundamental response of a bubble to a perturbation is 
a purely radial uniform volume pulsation in which the bubble 
oscillates about its equilibrium radius while maintaining its 
spherical shape. This approximates a perfect monopole source. 
Other responses, harmonics, may also occur, but the 
fundamental is the most efficient radiator of sound and 
dominates our interest. 
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The frequency of this fundamental volume pulsation is 
determined by the stiffness of the gas plus the surface 
tension and a mass contribution from entrained fluid. Many 
derivations have been published relating a bubble's 
equilibrium radius to its resonance frequency. One of the 
most complete treatments is found in Clay and Medwin [Ref. 3]. 
The resonance frequencies of air bubbles of given radii at 10 
cm depth of fresh water as calculated using Equation 1 [Ref. 
3:eq. 6.3.11:p. 197] are given in Table 1. Here, b and /3 are 
functions that correct for the adiabatic assumption and 
surface tension. These two corrections are counteracting and 
result in very little change to the value calculated by 
ignoring them. Nonetheless, they were included for accuracy. 



1 

f«= 

2na 



37bB Pa 2 



( 1 ) 



where, 

a = bubble radius in cm 

= ratio of specific heats, Cp/C^; taken to be 1.4 02 

= ambient pressure on the bubble corrected for 
depth in dynes/cm^ 

Pf^ = ambient density of the fluid surrounding the 
bubble taken to be 1.03 gm/cm^ 

See [Ref. 3] for expressions for b and f3 . 
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TABLE 1* RESONANCE FREQUENCY FOR A GIVEN BUBBLE RADIUS 
INCLUDING CORRECTIONS FOR DEPTH (10 CM), SURFACE 
TENSION AND THERMAL CONDUCTIVITY 



Bubble Radius, a fum) 

100 

70 

50 

40 

30 

20 

10 

9 

8 

7 

6 

5 



Resonance Frequency (kHz) 

31.4 

44.6 
62 . 1 
77.3 

102 . 6 
152.9 
303.8 

337.7 
380.4 

435.8 

510.8 
618.0 



* The data in this table is included in the Easy Reference 
Table in the Appendix. 
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Factors directly effecting the resonance frequency include 
the volume and type of gas, properties of the fluid medium and 
depth of the bubble. For this experiment the gas was 
generated by electrolysis. Various gases can be generated by 
this method depending on the concentration of chemicals, such 
as salts, in the water and the voltage applied. The gas 
produced in this experiment was probably hydrogen, Hj. 
Chlorine, Clj and Oxygen, Oj are the other possibilities. 
Since these gases are all diatomic, they share some common 
properties. The principle property of interest is the ratio 
of specific heats, which is listed in Table 2 for these gases 
and air. Also listed is the approximate error induced by 
using the specific heat ratio for air in each case. Since 
these errors are small and many bubble properties have been 
determined and published for air bubbles in water, air was 
assumed to be the gas in the bubbles for all calculations. 
This also allows these results to be more quickly related to 
any fresh water air bubble situation. 

TABLE 2 SPECIFIC HEAT RATIOS FOR SOME GASES 

GAS SPECIFIC HEAT RATIO RESONANCE FREQ. ERROR FROM AIR 

Air 1.402 0.00% 

Hj 1.41 0.28% 

Oj 1.40 0.07% 

CI 2 1.34 2.24% 
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other factors which effect the resonance frequency, such 
as the depth, pressure, radius and fluid density, are all 
interrelated. Assuming the mass of gas in the bubble stays 
constant, a depth increase raises the pressure and decreases 
the bubbles radius. All of these changes contribute to an 
increase in the resonance frequency. On the other hand a 
fluid density increase lowers the resonance frequency. A gas 
with a larger specific heat ratio would have a higher 
resonance frequency. For this experiment, all resonance 
frequencies are given for air bubbles at 10 cm depth of fresh 
water with a surface pressure of one atmosphere. 

The response of a bubble to an excitation is most dramatic 
when the frequency of the excitation is exactly equal to the 
resonance frequency of the bubble. The bubble's response at 
resonance is much greater than the response off resonance as 
shown in Figure 1 [Ref. 3;Figure 6.4.1: p.204]. 

Due to the large amplitude oscillations of the bubble 
driven at resonance, the higher order terms in the bubble's 
equation of motion are not small with respect to the linear 
terms. These higher order terms are necessary for the dual 
frequency method to work. 

B. DUAL FREQUENCY EXCITATION AND RESPONSE 

References 1 and 2 give extensive treatments of the 
oscillations of a bubble under the influence of two impinging 
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sound fields of different frequencies, f, and £ 2 - fi is taken 
to be a high frequency, much higher than the bubble's 
resonance frequency, yet with a corresponding wavelength 
longer than the diameters of the bubbles of interest. fj is 
a lower frequency, at or near the bubble's resonance 
frequency. As the bubble enters the sound fields, it is 
excited into large amplitude resonant oscillations by fj. It 
is also excited by f,, although at a much lower amplitude. 
The result is that the bubble extracts energy from these sound 
fields and reradiates sound as a monopole source. Because the 
large amplitude oscillations are nonlinear, the bubble 
radiates sound at f^ and f2, as well as at the sum (f^ + fj) 
and difference (f, - f2> frequencies. All other possible 

combinations are also radiated, but at much lower amplitudes. 
This natural frequency mixing is the cornerstone of the dual 
frequency method. 

1. The Resonant Bubble's Response to Dual Frequency 
Excitation 

Isolating our attention to the sum frequency, the 
pressure radiated by a bubble is given by Equation 2 [Ref. 1 ]. 

Pi P2 

P* = (2) 

p a 6 r 
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where. 



p = density of fluid 

= acoustic pressure amplitude at frequency f^ at 
the location of the bubble 

P 2 = acoustic pressure amplitude at frequency fj at 
the location of the bubble 

p^ = acoustic pressure amplitude at the sum frequency 
(f^ + fj) at the location of the receiver 

= resonance (angular) frequency of the bubble 
S = damping coefficient 
r = distance from bubble to receiver 

It has been assumed that fj is equal to f,,, the bubble's 
resonance frequency. The sum frequency pressure amplitude is 
inversely proportional to the distance, r, from the bubble. 
The damping coefficient, S, is claimed to be almost 
independent of frequency in Reference 1. Values of S for 
given bubble sizes at sea level are taken from Reference 3 and 
listed in Table 3. Using Tables 1 and 3, the density of the 
fluid and an assumed distance from the bubble to the receiver, 
the sum frequency pressure amplitude at the receiver can be 
calculated for the bubble sizes of interest in terms of p^ and 
Pj. This information is shown in Table 4. 

f^ performs a task similar to that of the carrier 
frequency of a radio transmission. Hence, in this thesis f^ 
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TABLE 3* DAMPING COEFFICIENT FOR SOME BUBBLE SIZES (NUMBERS 



READ FROM FIGURE 6.3.1, 
Bubble Radius, a (um) S 


P. 199, REFERENCE 3) 
(Damoina Coefficient) 


100 


0.075 


70 


0.085 


50 


0.095 


40 


0.100 


30 


0.115 


20 


0.13 


10 


0.16 


<10 


0.16 



* The data in this table is included in the Easy Reference 
Table in the Appendix. 
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TABLE 4* RATIO OF SUM FREQUENCY PRESSURE AT THE RECEIVER FACE 
TO THE PRODUCT OF SOUND FIELD PRESSURES AND P 2 AT 
THE BUBBLE 



p+ 

Bubble Radius, a fum) (x lO'^’Pa'’) 

P1P2 



100 

70 

50 

40 

30 

20 

10 

9 

8 

7 

6 

5 



43 . 2 

27.0 

17.4 

13.4 
8.79 
5.25 
2 . 16 
1.94 
1.72 
1.50 
1.28 
1.05 



* The data in this table is included in the Easy Reference 
Table in the Appendix. 
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is called the carrier frequency. It can be any convenient 
frequency high enough with respect to the resonance 
frequencies of the bubbles of interest that the harmonics of 
these bubbles are significantly less than the expected 
difference frequencies. This prevents confusing a harmonic 
with the difference frequency. The bubbles of interest must 
have diameters small with respect to the wavelength associated 
with f,, so that it oscillates uniformly. 

The bubbles of interest to this investigation have 
radii from 100 down to 7 ixm which corresponds to resonance 
frequencies from approximately 20 to 400 kHz. A reasonable 
minimum wavelength for these bubbles is 400 /xm, which 
corresponds to a maximum carrier frequency of 3.7 MHz. This 
ensures the acoustic wavelength is at least four times the 
largest bubble radius. Assuming harmonics higher than the 
third or fourth will be insignificant leads to a minimum 
carrier frequency of about 5 times the highest bubble 
resonance or 2.0 MHz. The actual value of f^ used in this 
experiment was 2.65 MHz which is near the center of this 
frequency range. 

The purpose of fj is to excite the bubbles into large 
amplitude resonant oscillations. For this reason, fj is 
called the excitation frequency. It must encompass the range 
of resonance frequencies of the bubbles of interest. This 
means that f 2 must cover the full range from 30 to 400 kHz 
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continuously. Two schemes for producing £2 are possible. 
One as sweeping the frequency repeatedly over the desired 
range. This method will periodically excite any bubble in the 
range, but only for a short time. Another option is 
broadcasting band limited white noise covering the proper 
frequency range. This more continuously excites any bubble 
in the range at its resonance frequency. This second method 
was chosen so that all bubbles will always find their 
resonance frequency in the fj sound field. This method 
requires more power since it broadcasts over a large bandwidth 
continuously at moderate amplitude. 

Using noise to ensure the resonance frequencies of 
interest are present also ensures that additional frequencies 
are present. This does not confound the dual frequency method 
though, since each bubble naturally responds strongest to its 
resonance frequency. Due to the large quality factors, the 
bubbles will vibrate strongly at their resonance frequencies 
and greatly ignore the sound energy at other nearby 
frequencies. This will allow the detected sum frequency to 
pinpoint the resonance frequency, f,,, of the bubble that 
generated it. f, = - f.grrier* 

As bubbles get smaller their acoustical cross sections 
(absorption, scattering and extinction) become smaller. To 
keep smaller bubbles excited strongly enough to produce 
observable dual frequency mixing, coloring the noise sound 
field is necessary. A ramp function that causes more sound 
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energy at the higher frequencies alleviates this problem 
somewhat. 

A relatively large bubble responds strongly at its 
resonance frequency and also scatters higher frequencies. The 
scattered amplitude can be higher than the amplitude due to 
a bubble resonant at the higher frequency due to the larger 
bubble's large scattering cross section. This leads to an 
ambiguity for other methods as shown by curve b of Figure 2 
[Ref. 1] . For a uniform distribution of bubble radii in the 
mixture, a high bubble resonance frequency may be detected at 
high amplitude due to either, (1) a resonant bubble or (2) a 
much larger scattering bubble. Bubble scattering detection 
techniques looking for the bubble's resonance frequency 
directly cannot be sure which of these mechanisms is causing 
the signal without more information. Reference 3 discusses 
this difficulty in more detail. 

The dual frequency method is not susceptible to this 
problem. Mixed frequencies, (f^ + f 2 ) and (f^ - f 2 ) , are 
produced only with the large amplitude whole body oscillations 
that occur around resonance. Therefore, since mixed 
frequencies are not simply a scattering phenomenon, it becomes 
unimportant that multiple scattering may be occurring. The 
mixed frequencies are dependent on the bubble's resonance 
frequency, f^, and become (f^ + f^) and (f^ - f^) . So a signal 
received at these frequencies indicates a bubble with 
resonance frequency, f„, is in the sound field. Even though 
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Figure 2 Comparison of the Dependence of Sum Frequency 
Pressure, p*, and Exciter Frequency Echo Pressure, 
Pj, on Bubble Radius, Curve a: p* and Curve b: p. 



the scattered signal amplitude peaks, dips and then rises in 
a ramp fashion above the peak as bubble radius increases, as 
shown in curve b of Figure 2, the sum frequency signal is 
strong for only one bubble size, as shown in curve a. 

2. Dual Frequency Interaction for Solids and Non-Resonant 
Bubbles 

Beyer [Ref. 4] has gathered a lot of information 
concerning the nonlinear interaction of dual frequency sound 
fields. In this compilation of results from his own work and 
others, he presents information on the generation of sum and 
difference frequencies by the interaction of the fluid medium 
alone, a solid in the fluid and a bubble in the medium. These 
treatments allow us to estimate a maximum expected signal for 
cases without resonant bubbles. 

a. Sound Interaction with Sound 

In an experiment by Jones and Beyer it was shown 
that no sum or difference frequency is generated when two 
sound beams cross at right angles. In fact, sum and 
difference frequency signals are generated only when the beams 
are collinear. This case is termed the "parametric end fired 
array." It results in a highly directional sound beam. The 
pressure of this sound beam is proportional to the pressures 
of the primary frequency signals and the frequency separation 
between these signals. To produce the signal, the sound beams 
must be collinear for some distance. As a rule this process 
is less than one percent efficient. 
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b. Sound Interaction with Solids 



The first effect to consider in the presence of 
a solid is scattering. The direction of scatter depends on 
the size, shape and orientation of the solid. To consider the 
parametric end fired array, assume that significant portions 
of both carrier and exciter are scattering toward the 
receiver. Although the receiver is very close to the sample 
volume, some interaction may occur during the interval. The 
resultant signal would have, at most, less than one percent 
of the scattered power. This power estimate is generous 
considering the short interaction length. 

The second effect of solids in a dual sound field 
is their motion due to the radiation forces of the two sound 
fields. The motion of the solid induces the sum and 
difference frequency in the medium it contacts. To have this 
effect, the sound fields must be significantly intense to move 
solids at both the carrier and exciter frequencies. 

c. Sound Interaction with Non-Resonant Bubbles 

Beyer gives theory and experimental results for 
an air bubble interacting with dual frequency sound. The 
frequencies used were so high (5 and 7 MHz) that the bubbles 
were not resonant at these frequencies. The theory considers 
the interaction of sound scattered by the bubble. Scattering 
is the primary interaction for bubbles larger than resonance 
size. For the case presented, a 35 /lim bubble is ensonified 
by 7 MHz at 2.45 x 10^ Pa and 5 MHz at 3.36 x 10^ Pa. These 
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frequencies are scattered omnidirectionally. The resulting 
pressures 48 cm away are 7 Pa and 6 Pa, respectively. The sum 
signal arrives with a pressure of 0.07 Pa. 

This is interesting compared to the same analysis 
for a solid sphere. Here the sum frequency pressure is 
proportional to the square of the solid particle's radius and 
varies with direction. The maximum sum frequency signal 
exists directly down range from the sphere on a line with 
either source transducer. For a sphere over 90 times larger 
than the bubble mentioned above and identically ensonified, 
the maximum sum frequency pressure was 0.007 P^ at these two 
points. The minimum signal was between these points at a 
pressure of 0.00007 Pa. The placement of the receiver between 
the two sources ensures that we avoid the maximum values of 
the scattered signal pressure. 

d. Overall Analysis of Solid Interference 

Clearly there is a limit to the amount of solid 
material that can pass through the dual frequency sample 
volume without great effect. Much, however, can be said for 
this method's insensitivity to solids. First, the above 
experiments were performed at greater drive pressures than 
required for the dual frequency method. Lowering the drive 
pressures lowers these nonlinear effects. Second, note that 
the solid has much less sum signal than the non-resonant 
bubble. The non-resonant bubble has much less sum signal than 
the resonant bubble, as shown in Figure 2. It follows that 
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a solid particle will have a sum signal that is much less than 
that for a resonant bubble. This provides the means for 
discrimination between bubbles and solids. 

The final argument to separate the resonant 
bubble's sum frequency signal from those of scattering 
phenomena is the frequency selectivity. The excitation signal 
is relatively low power over a broad band. Whereas the solid 
would scatter all frequencies uniformly, the bubble will be 
fundamentally excited by the presence of its resonance 
frequency. The whole body oscillations then mix the bubble's 
resonance frequency with the high power carrier. This sum and 
difference frequency signal will be easily distinguished from 
the other noise. 

C. ADVANTAGES AND DISADVANTAGES OP THE DUAL FREQUENCY METHOD 
The dual frequency method for bubble detection and 
resonance frequency determination has its good and bad points. 
The following is a quick summary of some of the key issues. 

* Advantages 

- conclusive resonance frequency determination 

- discrimination between solids and bubbles 

* Disadvantages 

complexity of implementation 

lower signal to noise levels than direct resonance 
frequency methods 
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* other Effects 



rectified diffusion 
radiation pressure 

1. Advantages 

The foremost advantage of dual frequency bubble size 
determination is the single peaked response shown in Figure 
2. Since this is not a simple scattered sound measurement 
technique, scatterers do not dominate the results. Large 
bubbles respond at their resonant frequency. They do not 
generate the same signals as small bubbles. This advantage 
is the same reason why solids do not give false detections. 

Since this is an acoustic method using actual 
resonance phenomena, true resonance frequencies of skin 
covered bubbles can be determined by experiment. This could 
be important for certain types of plankton or other resonant 
scatterers. The contribution of the skin need not be 
calculated or assumed. 

2. Disadvantages 

One disadvantage is the complexity of the experimental 
set up and support equipment. Function generators must be 
capable of producing colored noise and narrow band high 
frequency. Frequency spectral analyzers must be capable of 
quickly and repeatedly analyzing the received signal over a 
wide range of frequencies. The transducers must be placed 
correctly to overlap the sound fields and receive a signal. 
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The second major disadvantage is the relatively low 
signal to noise level. The amplitudes of the sum and 
difference frequency signals are very small compared to the 
amplitudes of the signals received at the carrier and 
excitation frequencies. The sum and difference frequency 
generation mechanisms are nonlinear interactions and weaker 
than the linear interactions. 

3. Other Effects 

Another problem might be that the sound field could 
change the nature of the bubble being sized. Rectified 

diffusion, a process by which an oscillating bubble actually 
grows, will change the resonance frequency. Then the detected 
resonance frequency would not be the original undisturbed 
frequency. Crum [Ref. 5] gives a detailed discussion of the 
subject of rectified diffusion. The pressure amplitude 

necessary to initiate rectified diffusion is [Ref. 5:eq. 19:p. 
217] 

(pa^Wn=)=[ (1-(W=/W«=) )=+6=W=/W„^](l+(2a/ap_)-Ci/c„] 

Ptm^— ( 3 ) 

(3 + 4K) (Ci./Co)-{[3(n-l) (3n-4)/4 J + (4-3n)K>[l+(2o/ap„) ] , 



where, 

w = frequency of excitation 
a = surface tension of the liquid 
p„ = ambient pressure 
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Cj = concentration of dissolved gas in the fluid far 
from the bubble 

Cq = equilibrium concentration of gas in the fluid 
K = thermal conductivity of the gas in the bubble 
n = polytropic exponent 

p^„ = threshold acoustic pressure at excitation 
frequency. 

Table 5 shows the results of some lengthy computations 
using Equation 3. Comparing the actual pressures used to 
these threshold values verified that rectified diffusion was 
not a problem with this experiment. The other saving point 
is that rectified diffusion takes hundreds of seconds to cause 
an appreciable change in bubble size as shown in [Ref. 5: Fig. 
7:p. 222]. This is much longer than bubbles are expected to 
remain in the sound field. 

Radiation force and streaming are two nonlinear 
effects discussed in Reference 4. These mechanisms have 
effects on the bubbles transiting the dual frequency sound 
field. The radiation force causes the bubbles to move away 
from the regions of high acoustic pressure. Streaming causes 
the bubbles to move away from the source transducers when in 
the sample volume. These motions make it difficult to 
maintain a small bubble in a small sample volume long enough 
to detect it. The effect of both of these mechanisms can be 
reduced by lowering the carrier and excitation pressure 
levels. Having a large sample volume with respect to the 
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TABLE 5. THRESHOLD PRESSURES FOR RECTIFIED DIFFUSION OF 
VARIOUS AIR BUBBLES IN AIR SATURATED WATER 

(q/C, = 1) 



Bubble Radius 
(/im) 


Resonance 
Frequency (kHz) 


Threshold 
Pressure (Pa) 


Frequency 

(kHz) 


100 


31.4 


1438 


31.4* 


50 


62 . 1 


2535 


62 . 1* 


10 


303.8 


8676 


303.8* 


5 


618.0 


1515000 


2650.0** 



* W/W^ = 1 

** = carrier frequency 
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distance a bubble moves under the influence of these forces 



would also mitigate their effect. 

D. BUBBLE RISE TIME 8I2IH6 

In order to verify that a bubble's acoustically indicated 
size is correct, some other method of determining the bubble 
size is necessary. A relatively easy, straightforward and 
accurate method is to measure the rise time of the bubbles as 
they rise through a known distance in still water under the 
influence of buoyancy. Reference 6 has documented this 
technique . 

Essentially, a bubble rises at constant speed due to the 
balance between its buoyancy and its drag. These forces have 
different dependencies on radius. Hence, the rise speed of 
bubbles varies with radius; the larger the bubble, the faster 
it rises. The terminal rise speed of a bubble is given by 
[Ref. 6] 

‘ (4) 

t 

where , 

U = rise speed 

g = acceleration due to gravity 
Cp = drag coefficient 



U = 



Sag 



3 Cl 
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The values determined for time to rise one inch for 
various bubble sizes are shown in Table 6. The Schiller and 
Nauman drag law was used for these computations as given by 
Equation 5 [Ref. 6:eq. 4:p. 6]. 

24 

Co = — (l+0.15R«°-«®^) (5) 



Rg = Reynold's number = 2aU/V^ 

Vf = kinematic shear viscosity of fluid 

In practice, there are a few difficulties with this 
method. The experimenter must be able to monitor the moving 
bubble. Thirty micron radius bubbles are about the smallest 
easily timed, with good lighting and clear water. There is 
some parallax since the scale cannot be extremely close to the 
rising bubble. With large bubbles, which move fast, observer 
response time is a factor. The experimenter cannot be 100 
percent certain that the bubble timed is not accompanied by 
other bubbles that effect the acoustic output. 

Overall, the rise time method is helpful. It 
allows at least a rough comparison of the acoustically 
determined resonance frequency with an expected value for 
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TABLE 6* TIME TO RISE ONE INCH FOR VARIOUS BUBBLE SIZES 



Bubble Radius, a ( um) 
100 
70 
50 
40 
30 
20 
10 
9 
8 
7 
5 



Time to Rise One Inch (sec^ 
1.44 
2.55 
4 . 61 
6.95 
12 . 0 
26.5 
105 
129 
163 
213 
416 



* The data in this table is included in the Easy Reference 
Table in the Appendix. 
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those bubbles that can be timed. Even when it cannot be used 
for every bubble, it gives the experimenter confidence in the 
dual frequency method. 
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PROBLEM APPROACH 



The goal of this thesis is to develop and test a system 
which uses the dual frequency method to detect and size 
bubbles. In order to accomplish this goal, the following 
capabilities had to be incorporated into the experimental 
system. 

* A method of generating bubbles and placing them in the 
sample volume. 

* A method of placing solid particles in the sample volume. 

* A method of producing a sound field of sufficient pressure 
amplitude. 

* A method of receiving the pressure signals from the bubble 
and determining the frequency spectrum of that signal. 

These methods are discussed in this chapter. 

A. BUBBLE GENERATION 

A method for generating bubbles of very small size is 
necessary for this study. Two basic generation methods were 
tried. In the first, bubbles were generated by pressurizing 
a pipette, made from a glass rod with a very narrow axial 
hole. The pipettes were very delicate, produced an audible 
sound during bubble generation and produced bubbles much 
larger than desired for this study. For these reasons, this 
method was abandoned and an electrolytic method was developed. 
Figure 3 shows the basic bubble production setup. 
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Tap 

Water 




Figure 3 Electrolytic Bubble Generation Setup 

Using the arbitrary (ARB) wave function mode on the 
Hewlett Packard (HP) 3314A, a negative electrical pulse was 
generated. The pulse was applied to a small wire, all but the 
very tip of which was insulated. A return wire with a bare 
coil tip completed the circuit back to the function generator. 
Ordinary (and by no means pure) tap water served as the 
electrolyte between the two leads. Very small bubbles would 
form and rise from the small wire during the negative pulses. 

The nature of the bubbles was easily varied, but not 
controllable. Because the small wire tip corroded and 
collected deposits rapidly with use, the bubble size and 
quantity would vary with time, even with all of the parameters 
held constant. The use of a non-corroding wire, platinum for 
instance, might have solved this difficulty. 

The number of possible variations offered by this setup 
proved useful. By adjusting the time scale on the ARB cycle 
the bubbles could be made to appear in quickly recurring 
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bursts having small numbers of small bubbles or in groups 
having a large number of larger bubbles with a larger interval 
between groups. As a rule to obtain smaller bubbles: (1) use 
a thinner wire, (2) use a shorter duration negative pulse, (3) 
use a lower voltage negative pulse. Conversely, larger 
bubbles were more likely with thicker wire, higher voltages 
and longer duration pulses. 

It was impossible to generate only one size bubble over 
any length of time with this method. With several small 
bubbles present, two would often join to form a larger bubble. 
The larger bubble would then rise faster, overtake and join 
another bubble. To prevent this, an on-off switch was placed 
in the circuit. The current was interrupted after a bubble 
was formed to limit interference with other bubbles. 

1. Variations of Bubble Control 
a. Grid 

Since the position of the rising bubbles is 
important to this experiment, some way is desired to move the 
generation point a known amount. One solution is a grid of 
32 wires terminating on a small blank circuit board. Arranged 
in 4 columns and 8 rows with the tip of each wire about one 
millimeter from its neighbor, the source of bubbles could be 
moved by switching one wire off and another on. Also a 
massive cloud of bubbles could be achieved by turning on all 
or some of the wires at the same time. 
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b. Bubble Size Discrimination Plate 

In order to narrow down the size range of bubbles 
which passed through the sound field, a small horizontal water 
jet was produced near the bubble generator. The jet was 
established by connecting a water reservoir to a small hollow 
glass tube. The speed of the jet could be controlled through 
the height of the reservoir. As bubbles were generated they 
were swept downstream by the jet. The length of travel was 
inversely proportional to the size of the bubble. Eventually 
the bubbles would rise out of the jet and continue up to the 
surface. A plastic plate with one small hole was placed above 
the jet. This arrangement allowed only those bubbles leaving 
the jet at the position of the hole to rise through the sound 
field, as shown in Figure 4. Although helpful in narrowing 
down the range of bubble sizes, some variation was still 
apparent . 




Wire 



To Water Tank 



To HP 3314A 



Figure 4 Bubble Size Discriminator Plate 
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c. Blow Pipe 

There are a number of difficulties associated with 
getting very small bubbles into the sound field. Worst of 
all, they are just about impossible to see. They drift away 
from the acoustic sample volume very easily, and they rise at 
very slow speeds, several minutes per inch. In order to get 
them into the sample space, they had to be convected. To this 
end, the bubbling wire was inserted into a narrow glass tube. 
The tube was pressurized by a water reservoir with about one 
foot of static head. The water was allowed to flow 
continuously while the bubbles came in short pulses. The tube 
could be (relatively easily) positioned so that the convected 
bubbles would pass through the sample volume. 

B. SOLID INJECTION 

In order to show that the frequency mixing due to solid 
particulates in the system is dramatically different from that 
due to bubbles, and that therefore the system would show the 
ability to discriminate between a bubble and a solid, such 
solid particles had to be introduced into the sample volume. 

Three basic types of solid particles were injected. The 
first, dry sand, was dropped directly over the dual sound 
field via a dry funnel and allowed to fall to the tank bottom. 
It was never certain whether the particles fell directly 
through the most intense portion of the sound field or not, 
due to their long fall and fluttering motion. 
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Wet sand was injected in a water slurry via a long wet 
funnel terminating just above the dual sound field. Wetting 
the sand prior to injection reduced the number of visible 
bubbles entrained with the sand. 

Thin wall hollow glass beads were also injected into the 
sound field. Due to their extreme buoyancy they had to be 
convected. A glass tube was aimed at the dual sound field 
horizontally and pressurized with a bead slurry. The beads 
would pass into the sound field horizontally and rise to the 
surface. The glass beads ranged in size from about 5 to 50 
Mm radius, as determined by microphotography. The sand 
particles ranged from 1 to 1 mm in largest dimension. 

C. TRANSDUCER CHARACTERISTICS AND CONFIGURATION 

Due to the range of frequencies used in this experiment 
several types of transducers were needed. The lower frequency 
driving transducer, called the exciter, needed significantly 
different parameters than the higher frequency driving 
transducer, called the carrier. Since the sum and difference 
frequencies were close to the carrier frequency, the receiving 
transducer, or receiver, was very similar to the carrier 
transducer. Other transducers were used for special purposes, 
such as calibration, but were not part of the dual frequency 
system. This section will present each transducer's 
characteristics and the methods used to determine those 
characteristics . 
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1. Types of Transducers 



Two 6 cm diameter focused piezoelectric transducers 
were used for the carrier and receiver. They were chosen for 
the acoustic gain advantage of focusing. The acoustic 
pressure at the focus is many times that at the face of the 
transducer. The serial numbers of these transducers are C5575 
and C5574. They were manufactured by Harisonic. Nominal 
specifications include a 2.25 MHz resonance frequency and a 
3 inch focal length. 

A Panametrics V301 3 cm diameter circular 
piezoelectric disk transducer was used as the exciter. This 
choice was made since it had good frequency response over the 
range of 30 to 400 kHz. This transducer had a source strength 
that mainly increased with frequency over this range. That 
was ideal for this purpose. The serial number of this 
transducer was 93598. To assist in reciprocity calibration 
a similar transducer, serial number 93596, was used. 

To determine the transducer beam patterns and assist 
in correctly positioning the transducers in the system, a very 
small hydrophone, made by Specialty Engineering, was used. 
The hydrophone looks like the tip of a very sharp pencil, with 
all dimensions of the ceramic receiver less than a millimeter. 
A preamplifier supplies enough signal boost to give 
significant oscilloscope readings. This hydrophone was used 
to probe the sound fields for relative measurements. No 
attempt was made to calibrate the probe. 
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Finally, an EDO corporation model 6600 spherical 
piezoelectric hydrophone was used for the reciprocity 
calibration of the exciter. 

2 . Beam Patterns 

The beam patterns of the exciter, carrier and receiver 
are important since their intersection determines the size of 
the sample volume. Also, an accurate knowledge of the beam 
pattern is crucial to the calibration of the focused 
transducers. The methods used to measure the beam patterns 
of the transducers are described in this section, 
a. Beam Pattern Method 

In order to obtain an accurate beam pattern, the 
transducer and probe hydrophone were mounted to a common frame 
via a network of micrometer positioners. Thus, the probe 
could be moved known distances, relative to the transducer, 
in small steps. The setup allowed relative motion in the 
three orthogonal directions of vertical, range and cross 
range. 

The first step in measuring the beam patterns of 
the focused transducers was to locate the position of the 
focus. Then readings of distance and probe output voltage 
were taken relative to the focus. At any given frequency, the 
probe output voltage was assumed to be proportional to the 
acoustic pressure at the probe. Normalized pressure values 
were obtained by dividing each value in a data set by the 
maximum value in that data set. 
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The exciter beam pattern was much broader than 
those of the focused transducers. Cross range data for the 
exciter was taken at a few ranges of interest. Since the 
exciter's beam pattern is, theoretically, narrower at higher 
frequencies, the data was taken at 400 kHz. This verified 
proper coverage at the highest exciter frequency of interest. 
The patterns at lower frequencies were then guaranteed to be 
sufficiently broad. 

b. Beam Pattern Results 

At 400 kHz the exciter had a very smooth beam 
pattern. At a distance of 4.0 cm the 3 dB beam width was 
about 1.0 cm. Figure 5 shows the normalized pressure versus 
cross range at a range of 4.0 cm. The variation of pressure 
with range was minimal for ranges 4 cm to 6 cm. At a range 
of 10 cm the maximum pressure was approximately 3 dB less than 
that at 4.0 cm range. Figure 6 shows the area within which 
the pressure was greater than the 3 dB less than the maximum 
value (0.71 Pn«x) • As will be shown below, this beam pattern 
is sufficient to ensonify the intersection volume of the 
carrier and receiver, which limited the size of the sample 
volume. 

The beam patterns of the focused transducers were 
much more complicated. The acoustic pressure changed rapidly 
with distance, especially at the focus. Figure 7 shows the 
variation of normalized pressure with range along the axis for 
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Nomnallzed Pressure 




Figure 5 Exciter Cross Range Pressure 
Variation at a Range of 4.0 cm 
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Cross Range (mm) 




Figure 6 View of Exciter Main Beam. Area Inside curves Is 
High Pressure Region. 
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Normalized Pressure 




Figure 7 Focused Transducer's Pattern of Pressure versus 
Range at the Focus 
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transducer C5574. 



Both focused transducers were similar in 



this regard. At the focus, the pressure stays at or near the 
maximum for about 5 mm of range, and at least half of the peak 
value for 10 mm. Figures 8 and 9 show the vertical variations 
of normalized acoustic pressure at the focus of C5574 and 
C5575, respectively. The pressure stays at or near its 
maximum value for only 1 mm vertically. Comparing these 
graphs for the two transducers reveals that C5575 focuses a 
little tighter than C5574. The cross range beam patterns are 
very similar to the vertical patterns and are not shown. 
Although the beam patterns are axially symmetric in general, 
there are a number of small variations in the pressure field. 

Figures 10 and 11 show the horizontal plane (cross 
range versus range) containing the focal axis for C5574 and 
C5575, respectively. Isobaric lines surround areas within 
which the pressure is at least as great as the isobaric 
pressure. This is just another way to show that pressure 
varies ten times faster with cross range than with range. To 
compare the graphs, note that although the isobaric lines are 
labeled differently in units of dB relative the maximum value 
for that plot, the inner oval for C5574 represents the same 
pressure as the inner oval for C5575. These plots were 
instrumental in understanding the focused transducers' sound 
fields. They provided key information for the calibration 
processes as shown in the next section. Further detailed 
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Normalized Pressure 




Figure 8 C5574 Vertical Variation of Pressure at the Focus 



41 



Normalized Pressure 




Figure 9 C5575 Vertical Variation of Pressure at the Focus 
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Cross Range (mm) 




Figure 10 C5574 Horizontal Focal Plane Pressure Distribution 
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Cross Range (mm) 




Figure 1 1 C5575 Horizontal Focal Plane Pressure Distribution 
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analysis of the resulting sample volumes is given after the 
calibration section. 

3. Transducer Calibration 

The beam patterns in the previous section were all 
plotted relative to some maximum value for that transducer. 
To apply the dual frequency theory and calculate the expected 
sum and difference frequency signal pressures, the carrier and 
exciter pressures must be determined. Also the receiver 
sensitivity must be known. To this end, several calibration 
methods were used. This section presents a brief outline of 
each method and the results obtained, 
a. Methods of Calibration 

(1) Reciprocity. References 7 and 8 present a 
significant amount of information on reciprocity calibration 
of transducers. If done correctly, this type of calibration 
is a prime standard of transducer performance. The key to the 
reciprocity method is two fold, as Reference 9 points out. 
First, for a linear, passive and reciprocal transducer there 
is a simple relationship between the transducer's receiving 
sensitivity, M, and its transmitting response or source 
strength, S; 



M 

- = J, 

S 



( 6 ) 
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where , 

M has units of voltage/pressure 
S has units of pressure/current 
J has units of (voltage) ‘ (current) / (pressure) 

J is called the reciprocity parameter and depends on the 
nature of the calibration method. Several methods of 
reciprocity calibration are presented in Reference 8 with 
their associated reciprocity parameter. The second key to the 
reciprocity method is that the beam pattern of a transducer 
is the same whether it is acting as a source or receiver. 
This property was assumed in order to determine the beam 
pattern of the receiver in the last section. 

The full reciprocity method requires one 
transmitter, T, one hydrophone, H, and one reversible 
transducer, R. Three basic steps exist in the experiment. 
First, place the transmitter, T, a known distance, d, from H. 
Record the open circuit voltage response from H, Vht, for a 
known current, it into T Second, replace H with R. Record the 
open circuit voltage out of R, Vrt, for the same current into 
T. Third, replace T with R and return H to its original 
position. Record the open circuit voltage out of H, V„r, for 
a known current, 1 r, into R The hydrophone sensitivity, Mr, 
can be calculated according to 
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1 

2 



( 7 ) 






Vht Vhr 

J 

Vrt 1r 



where, 



J = 2\d/p^c, for spherical spreading 

\ = acoustic wavelength 

= characteristic impedance of the medium. 

A transducer's source strength, S, is a measure of how much 
pressure, p, it produces at a distance of one meter for a 
given current, i. When p is measured at a distance, d, 

P d 

S = (8) 

i 1 (m) . 



Sensitivity, M, is a measure of the open circuit voltage, V, 
produced for a given pressure, p, at the transducer face, 

V 

M = - (9) 

P- 

For this reciprocity calibration it is assumed that the 
hydrophone is small enough to not disturb the sound pressure 
field. Thus, the pressure in the above two equations can be 
equated, giving 
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s 



V d 



( 10 ) 



Mi 1 (m) . 

This equation is used to calculate the exciter's source 
strength after the EDO hydrophone's sensitivity is determined. 

(2) Two Transducer Reciprocity. A short cut can 
be made if two identical reciprocal transducers are to be 
calibrated. Here, only one transducer arrangement is 
required. The transducers are placed a known distance, d, 
apart. Then the open circuit voltage of one is measured while 
a known current is put through the other. The sensitivity and 
source strength of both transducers are given by Equations 11 
and 12 . 

M 

and S = - (11 and 12) 

J 

The disadvantage of this method is that the transducers must 
be exactly identical in sensitivity and source strength. This 
must be shown by another method before the two transducer 
reciprocity becomes valid. 

(3) Self Reciprocity Method. A third variation 
on the reciprocity method involves only one transducer. This 
method is very similar to two transducer reciprocity, except 
that there is no longer the need for identical transducers. 
The transducer is aimed at a perfect reflector, such as an 
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air-water interface. A pulse of known current is put through 
the transducer. The acoustic pulse is reflected from the 
interface and returns to the transducer. The open circuit 
voltage amplitude of the received pulse is recorded. Care 
must be taken to ensure the transducer is aimed correctly to 
receive the reflected pulse. The source strength and 
sensitivity equations are identical to those for the two 
transducer methods. Since the transducer remains connected 
to the pulse generator while measuring the return voltage, a 
diode or other high impedance arrangement must be used to 
isolate the influence of the pulse generator's low impedance 
on the open circuit voltage measurement. 

(4) Radiation Force Target Deflection. In this 
method, a perfectly reflecting target is hung in the 
transducer's sound field. The target is deflected due to the 
radiation force. The amount of deflection is determined by 
the' balance of the radiation force and the weight of the 
target. As shown in Reference 10, the total acoustic power, 
n, acting on the target is determined by the deflection 
distance, e, multiplied by a target specific factor, as shown 
in Equation 13. 



mgc 

n = e 

L 



(13) 
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where , 



m = mass of the target 
g = gravitational acceleration 
c = speed of sound 
L = length of target suspension 
e = deflection distance. 

At the focus, the total acoustic power is given 
by Equation 14. 



n 



I F dA F 



(14) 



where , 

Ip = intensity of the sound over the infinitesimal 
area, dAp 

dAp = infinitesimal area in a plane intersecting the 
focus and perpendicular to the transducer's axis 

The intensity is related to pressure, for plane or spherical 
waves, by Equation 15. 



so. 



I = (15) 

PoC 



IT = 



J oC 



dAr 



(16) 
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Now, define s(x) as the pressure's radial 
shape function which equals the normalized pressure at a 
vertical or cross range distance, x, from the focus. In other 
words, this shape function looks exactly like Figure 8 for 
C5574 or Figure 9 for C5575. The pressure at a radial 
distance, x, in this focal plane, Pp(x), can be written in 
terms of the maximum pressure and the radial shape function 
as in Equation 17. 

Pf(X) = PFMA>tS(X) (17) 



substitution gives. 



n= 2n s^(x) xdx (18) 

PoC Jo 

PrM/vx ^ 

n = (19) 

PoC 

Agff is the effective beam area at the focus. It is obtained 
by numerically integrating the shape function squared over the 
focal plane. 

Substituting the target deflection parameters 
for n and solving for maximum focal pressure yields Equation 
20 . 






PoC^ mge | 



Aa-e -e L 



( 20 ) 



51 



In order to compare this result to that of 
the self reciprocity method, the pressure at the transducer 
face (determined from the source strength) must be converted 
to the pressure at the focus. The conversion is made assuming 
the acoustic energy at the face equals that at the focus. 
This assumption yields 



where , 



Pfmax 



P T ^ At j 

Ae« * 



= area of the transducer face 
Pj = pressure at the transducer face. 



( 21 ) 



b. Calibration Results 

Some of the methods described in the previous 
section were more suitable than others depending on the 
particular transducer to be calibrated. For instance, the 
full reciprocity method was ill suited for the focused 
transducers because of the rapidly varying pressure pattern 
in the focused sound field and the high frequency. Under 
these circumstance, the assumption that the hydrophone does 
not disturb the pressure field is not valid. The other tree 
methods were viable alternatives. For the exciter, the Self 
Reciprocity and Radiation Force Target Deflection methods were 
not used because of the low pressure amplitudes. 
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(1) Exciter Calibration. A full reciprocity 
calibration was performed on the exciter, using a second 
transducer of the same design (twin) and an EDO hydrophone. 
This type of calibration results in a value for the 
hydrophone's receiving sensitivity. The exciter's source 
strength can then be calculated by Equation 22. 

VedO d 

Se = ^ (22) 

Medo ie l(m) 



where, 

Sg = exciter's source strength 
Vgpo = open circuit voltage of EDO hydrophone 
MgpQ = EDO hydrophone's receiving sensitivity 
ig = current input to exciter 
d = distance between source and hydrophone 

The calibration was performed twice at each 
frequency to get two values for the EDO's sensitivity. To 
make the two values somewhat independent the roles of 
projector and reciprocal transducer were exchanged between the 
exciter and its twin. The values for the EDO's sensitivity 
were averaged and compared with expected values provided from 
the manufacturer. The average EDO sensitivity was used in all 
computations . 

The source strengths of the exciter and its 
twin were calculated using Equation 22. At some frequencies 
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the source strengths of the exciter and its twin were close 
to identical. The two transducer reciprocity calibration was 
performed twice in that case. Again, the roles of source and 
receiver were switched to provide some small degree of 
independence between the two runs. For this method the source 
strength is given by Equation 23. 

1 

2 



V = open circuit voltage of the receiver 
i = input current to the source 
J = reciprocity parameter, 2ld/p^c 
X = acoustic wavelength 

For these calibrations, the received open 
circuit voltage was measured directly on an oscilloscope 
having a 1 megaohm input impedance. The input current was 
determined by measuring the voltage across a 10.1 ohm 
resistor, placed in series with the driving transducer, on 
another channel of the same oscilloscope. The drive signal 
was obtained from a Hewlett Packard (HP) 3314A Function 
Generator. 

The full reciprocity method seemed to be the 
most reliable of these calibrations. The full reciprocity 
calibration results for the exciter were used in all further 



S = 



V 



i J 
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equations requiring exciter source strength. Since the 
exciter was used for the full continuum of frequencies from 
30 to 400 YB.Z , it was calibrated at 50 kHz intervals from 50 
to 400 )cHz. Table 7 contains the results of these 
calibrations . 

The true matter of interest is the pressure 
in the sample volume, which is obtained from Equation 8 
rearranged as p = s i'l(m)/r. The drive current is determined 
by measuring the voltage across the 10.1 ohm resistor .as 
discussed earlier. However, the drive signal is broadband 
noise. The oscilloscope cannot be used to determine voltage 
at a chosen frequency in this case. The voltage (rms) was 
measured with the HP 3585A Spectrum Analyzer. The peak 
current is therefore given by Equation 24. 

v/2 

i = (24) 

R 

Using this expression the pressure is found to be, 



P 



■J2 Vrmo S‘l(m) 
R r 



(25) 
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TABLE 


7 EXCITER CALIBRATION RESULTS. 

SOURCE STRENGTH IN UNITS OF Pa/mA 

FULL RECIPROCITY TWO TRANSDUCER METHOD 


1 . fkHz^ 


SfEXCITERI 


SfTWIN) 


S fEXCITERI 


S fTWINI 


50 


310 


526 


— 


— 


100 


1440 


1160 


-- 


— 


150 


1860 


2510 


— 


— 


200 


3340 


3250 


3380 


3440 


250 


5160 


5300 


5400 


5460 


300 


11100 


5560 


— 


— 


350 


11300 


10300 


11200 


11200 


400 


13600 


13700 


14200 


14200 
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where, 



= voltage determined by spectrum analyzer 
R = resistance of series resistor 

r = distance from the transducer face to the sample 
volume 

Table 8 gives the values of peak pressure for 
the maximum noise signal used to drive the exciter during the 
experiments. It turns out that the exciter has some 
significant variation in output pressure over this frequency 
range. Interpolated values were used for frequencies of 
interest that were not part of the calibration. 

(2) Focused Transducer Calibration. The first 
method attempted in the focused transducer calibration was the 
two transducer method. The position of both transmit and 
receive focused transducers were adjusted to maximize the 
return signal. This position turned out to be such that the 
transducers faced each other with the acoustic foci 
overlapping at half the distance between them, as shown in 
Figure 12. The two transducer method was performed twice. 
Switching the roles of source and receiver between C5575 and 
C5574. The obvious problem with this method was that it gave 
no indication of differences between the transducers. The 
result is a single number for source strength and another for 
sensitivity. It was suspected and later proven that the 
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TABLE 8* PEAK EXCITER PRESSURE AT 4 CM AT VARIOUS 
FREQUENCIES WHEN DRIVEN WITH MAXIMUM NOISE 
SIGNAL 



Frequency fkHz^ Peak Pressure at 4 cm (Pa^ 



50 


115 


100 


261 


150 


244 


200 


370 


250 


515 


300 


980 


350 


1000 


400 


1190 



* The data in this table is included in the Easy Reference 
Table in the Appendix. 
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Figure 12 Two Transducer Reciprocity 
Focused Transducers 



Calibration 



Using 
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transducers were not identical. This calibration resulted in 
a ball park verification of other methods used. 

Since the impedance of the focused 
transducers was not large compared to 10.1 ohms, and external 
resistor was not used in the determination of input current. 
Instead the magnitude of the free field impedance, |z|, of the 
source transducer was measured separately at the frequencies 
of interest using a HP 4192 LF Impedance Analyzer. With this 
value known, the input current could be deduced from the drive 
voltage. The receiving transducer's open circuit voltage was 
measured directly on the oscilloscope. 

As stated previously, one key to any 
reciprocity calibration is the reciprocity factor, J. For 
small, omnidirectional sources in open water this factor is 
normally given by Equation 26. 



2Xd 

J = (26) 

where, 

A = acoustic wavelength 
PoC = characteristic impedance of the medium 

The distance factor, d, is included to 
account for losses due to spherical spreading. However, for 
the focused transducers, spherical spreading is not the case. 
The sound field converges to and diverges from a focal point. 
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On the average, this can be equated to a plane wave phenomena 
in the sense that the receiver intercepts essentially all of 
the acoustic energy produced by the transmitter. Reference 
8 gives a brief presentation of plane wave reciprocity. 
Reference 11 provided a verification of the plane wave 
reciprocity factor given in Reference 8 and Equation 27. 

2Xd 

= (27) 

P^c 



where , 

A = area of transducer face (curved area for focused 
transducers) 

In order to avoid the need for identical 
transducers, the self reciprocity method was used. The 
transducer was positioned vertically underwater to achieve a 
maximum surface reflected signal. This placed the focal point 
at the air water interface as shown in Figure 13. 

A pulse of about 200 cycles at the desired 
frequency in the range of 2.25 MHz to 3.10 MHz was generated 
by the HP3314A every 10 milliseconds. The drive pulse was 
powerful enough to pass through the diode network and drive 
the transducer under test. However, the diodes caused some 
minor signal distortion. To account for this, the drive 
signal was monitored on a spectrum analyzer at each drive 
frequency used. A ratio of signal voltage at the drive 
frequency as measured on the spectrum analyzer to the waveform 
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Figure 13 Self Reciprocity Method for Focused Transducers 
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peak to peak voltage as measured on the oscilloscope was 
calculated for each drive frequency. This ratio was used as 
a correction factor to the drive voltage measured on the 
oscilloscope. 

The diodes essentially eliminated the HP 
3314A's low impedance from the receiving circuitry. Because 
of losses to the diode circuit, the received voltage measured 
on the oscilloscope with the diode circuit installed was 
noticed to be about 99 percent of the open circuit voltage. 
A factor of 1.01 was used to correct the oscilloscope readings 
to open circuit voltage readings. Since the received signal 
was spectrally dominated by the desired frequency, a second 
frequency correction factor was not required. 

Table 9 shows the results of this 
calibration. In general, C5575 was a more powerful source and 
a more sensitive receiver than C5574. Notice that the two 
transducer method result of 0.135 mV/Pa sensitivity at 2.65 
MHz is not far from an average of the self reciprocity method 
results for the same frequency. 

Note also the unexpected difference in the 
characteristic rolloff from the nominal resonance frequency 
of 2.25 MHz. C5575 shows an expected pattern of falling 
sensitivity with increasing frequency beyond this point. 
While C5574 stays virtually flat over the sum and difference 
frequency range. For this reason, C5574 was chosen to be the 
receiver. C5575 provided the carrier signal. 
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TABLE 9 


SELF RECIPROCITY CALIBRATION RESULTS FOR THE 
FOCUSED TRANSDUCERS 






C5575 




C5574 


FREO. 

(MHz) 


** SOURCE 
SENSITIVITY STRENGTH 

(mV/ Pa) (Pa/mA) 


♦SENSITIVITY 
(mV/ Pa) 


♦♦SOURCE 

STRENGTH 

(Pa/mA) 


2.25 


0.244 


58.3 




0.0982 


23 . 5 


2.30 


0.238 


57.0 




0.0969 


23.2 


2 .35 


0.236 


56.6 




0.0955 


22.8 


2.40 


0.232 


55.5 




0.0951 


22.7 


2.45 


0.232 


55.5 




0.0939 


22 . 5 


2.50 


0.230 


55.0 




0.0936 


22.4 


2.55 


0.230 


55.1 




0.0932 


22.3 


2.60 


0.229 


54 . 8 




0.0934 


22 . 3 


2 . 65 


0.227 


54.4 




0.0932 


22.3 


2.70 


0.226 


54 . 1 




0.0937 


22.4 


2.75 


0.224 


53.6 




0.0948 


22.7 


2.80 


0.221 


52 . 9 




0.0947 


22 . 6 


2.85 


0.217 


52.0 




0.0967 


23 . 1 


2.90 


0.213 


50.9 




0.0966 


23 . 1 


2 .95 


0.207 


49 . 4 




0.0982 


23 . 5 


3 . 00 


0.196 


46.9 




0.0974 


23 . 3 


3 . 05 


0.186 


44.5 




0.0972 


23.2 


3 . 10 


0.176 


42.1 




0.0966 


23.1 


* This 


data is 


included in the 


Easy 


Reference 


Table in the 



Appendix. 

** at the transducer face 
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of particular interest is the pressure of the 
carrier signal in the sample space. The maximum pressure at 
the focus, Pfhax/ is calculated from Equation 28 after 
determining the input current, i, and the effective acoustic 
area, A^ff. 






(S2.6smm== X i)= A-r 






(28) 



i 



V 






2 I Z I a 



e 5MHZ 



(29) 



where, 

Vpp = pea)c to pea)c voltage of drive wave form on the 
oscilloscope 

|zl2 65HHz magnitude of free field impedance of C5575 at 2.65 
MHz 

S 2 . 65 HHZ = source strength of C5575 at 2.65 MHz 

Since the carrier is monofrequency, the only 
parameter that varies is Vpp. Table 10 shows some calculations 
for PfHAx based on some typical input voltages. The following 
values were used for this table. S 2 . 65 hhz = 54. Pa/mA, A-, = 
3094.5 mm , Aeff = 0.6811 mm , |z|2,6shhz = 135.3 ohms. 
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TABLE 10* MAXIMUM PEAK FOCAL PRESSURE OF THE CARRIER AT 2.65 
MHz FOR VARIOUS INPUT VOLTAGES AS DETERMINED BY 
THE SELF RECIPROCITY CALIBRATION 



Voltaae fPeak to Peak) 


Maximum Peak Focal Pressure fPa) 


0.052 


707 


0.165 


2242 


0.41 


5570 


0.55 


7472 


1.10 


14940 


2 . 08 


28260 


2.55 


34640 


4 . 05 


55020 


5.2 


70650 


8.2 


111400 


10.35 


140600 


12.9 


175300 



* The data in this table is included in the Easy Reference 
Table in the Appendix. 
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Figure 14 Radiation Force Target Deflection Calibration Setup 
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An additional calibration using the radiation 
force target deflection method was performed in order to 
verify the reciprocity calibration results. This test was 
easily performed using the setup shown in Figure 14. The 
target was positioned at the approximate location of the focus 
in order to intercept the entire signal. The deflection was 
measured with a telescope mounted to a micrometer positioner. 
The HP 8640B Signal Generator created the 2.65 MHz signal with 
very fine voltage control and spectral purity. Using Equation 
20, Table 11 was generated using the deflection distances and 
input voltage data from this experiment. This determination 
of peak carrier pressure at the focus compares favorably with 
the self reciprocity method. 

4. Test Setup and Sample Volume 

Due to the highly focused beam patterns, precise 
alignment of the transducers was required. For this reason, 
the transducers were mounted to micrometer positioners which 
allowed adjustments of range, cross range and height. The 
alignment procedure consisted of first placing the probe 
hydrophone in the focal region of the carrier. The receiver 
was then carefully positioned so that its focus also coincided 
with the position of the probe. Similarly, the exciter was 
aimed at the probe to provide maximum excitation signal to the 
sample volume. The bubble generator was then positioned such 
that the bubbles impinged on the probe and thus would go 
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TABLE 11* 



MAXIMUM PEAK FOCAL PRESSURE OF THE CARRIER AT 2.65 
MHz FOR VARIOUS INPUT VOLTAGES AS DETERMINED BY 
THE RADIATION FORCE TARGET DEFLECTION CALIBRATION 



Voltaaes fPeak to Peak) 


Maximum Peak Focal Pressure fPa) 


4.0 


47080 


5.0 


66580 


6.0 


81540 


7.0 


88080 


8.0 


105300 


9.0 


115300 


10.0 


133200 


11.0 


148900 


12.0 


166400 


13.0 


188300 


14.0 


199700 


15.0 


215700 


16.0 


230600 


20.0 


264200 


25.0 


339500 


30.0 


410400 



* The data in this table is included in the Easy Reference 
Table in the Appendix. 
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through the sample volume. The probe was then removed and the 
experiment proceeded. 

Figure 15 is a full size drawing of a focused 
transducer and its focal region. This representation 

emphasizes how small the focal region is for the carrier and 
receiver. Figure 16 is a similar figure for the exciter. 
Figure 17 shows the three beam patterns overlayed (not to 
scale) . The intersection of these three regions is the sample 
volume. The receiver and carrier can be aligned no closer 
than 55 degrees due to the size of the transducers. The 
exciter provides full coverage when properly aimed from almost 
any angle. A value for the sample volume size can be 
determined for a given minimum acoustic pressure assuming the 
receiver and carrier foci overlap. The region can be 

approximately pictured in three dimensions as the intersection 
of two circular cylinders at an angle of 55 degrees. The 
radius of the cylinders depends on the acoustic pressure 
required for the dual frequency method to be effective. In 
the horizontal plane containing the focus, the area, 
through which a bubble must rise to be in the sample volume 
is given by 



4 To Tr 

= (30) 

sin 0 
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Figure 16 side view of the Exciter Transducer Showing the 
High Intensity Portion of the Beam Pattern (Actual 
Size) 





Figure 17 Beam Pattern Intersection 



where, 



= effective radius of carrier cylinder 
r^ = effective radius of receiver cylinder 
9 = beam pattern intersection angle 

The maximum vertical dimension of the sample volume is the 
minimum of the diameters of the intersecting cylinders, called 

^MAX • 

These equations were used to calculate the values 
shown in Table 12. The sample volume is not isobaric. The 
pressure is maximum in the center. Volumes of decreasing 
pressure surround the focus much like the skin of an onion. 
The geometry gets much more complicated farther axially from 
the focus, especially if the foci are not at the same point. 

Using the known sizes and rise velocities of the 
bubbles, it is possible to determine the length of time a 
bubble is in the carrier field at a given pressure. It is 
also important to realize that it is very unlikely that the 
bubble will pass directly through the focus. Figure 18 is a 
simplified view of the sample volume cut into volumes of 
relatively constant pressure. Top and side views are shown. 
The side view is looking axially toward the carrier. The 
three dotted lines show single paths by which bubbles may rise 
through the sample volume. The paths are labeled by the 
maximum carrier signal pressure encountered. The exciter 
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TABLE 12 SAMPLE VOLUME SIZE CHARACTERISTICS FOR 0 = 55° 



CARRIER 



PRESSURE 


r = r 

c r 


^AX 


AA 


^MAX 


ah 


(dB 


re 


MAX) 


(im) 


(Inin‘S) 


(rilin'^) 


(rain) 


(mm) 


0 


to 


-1 


0.2 


0.20 


0.20 


0.4 


0.4 


-1 


to 


-3 


0.3 


0.44 


0.24 


0.6 


0.2 


-3 


to 


-6 


0.5 


1.22 


0.78 


1.0 


0.4 


-6 


to 


-14 


0.8 


3.13 


1.91 


1.6 


0 . 6 


-14 


to 


-20 


1.0 


4 . 88 


1.75 


2 . 0 


0.4 


-20 


to 


-35 


1.5 


10.99 


6.11 


3 . 0 


1.0 
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Figure 18 Simplified View of the Sample Volume 



field is assumed constant and maximum over this small area. 
Table 13 depicts the time history of 100 lira, 20nm, and 9 iJ.m 
bubbles as they rise along these three paths. It is clear 
that a 30 kHz bubble stays in the sample volume no longer than 
0.17 seconds. Though this time is enough for the bubble to 
be excited it is very short for the receiving equipment to 
detect the signal. If rise time was the only consideration, 
then any bubble smaller than 30 /im radius would be in the 
sample volume at least 1 second. Unfortunately, radiation 
force of the exciter and carrier push bubbles out of the 
sample volume. The smaller the bubble, the more it is pushed 
and less likely to stay in the sample volume. To reduce this 
effect, the maximum pressure must be reduced. This, in turn, 
reduces the sample volume further. A compromise must be made. 

D. ELECTRONIC EQUIPMENT 

This section will detail the primary experimental setup 
in two major sections, drive signal generation and receiving 
equipment setup. The total system diagram is shown in Figure 
19. 

1. Drive Signal Generation 

The carrier signal was a 2.65 MHz sine wave. Because 
the sum and difference signals flank the carrier on a spectral 
plot, any anomalies near the carrier frequency could not be 
tolerated. This criterion led to the selection of the HP 
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TABLE 13 TRAVEL TIME AND CARRIER PRESSURE FOR BUBBLES 
TRAVELING ALONG THE PATHS SHOWN IN FIGURE 18 



TIME (mSEC) 


PRESSURE 


(dB re 


MAX) 


BUBBLE 100 ixm 20 


fim 9 )um PATH OdB 


-6 dB 


-20 dB 


0 


0 


0 


-35t 


— 


— 


5 


89 


434 




-35T 


— 


22 


398 


1939 






-35t 


28 


521 


2538 


-20t 






31 


660 


3216 




-20t 




40 


729 


3552 


-14t 






50 


912 


4444 




-14t 




57 


1042 


5077 


- 6T 






68 


1250 


6091 


- 3t 






74 


1354 


6598 


- It 






85 


1563 


7616 


0* 


- 6* 


-20* 


96 


1771 


8629 


- 14 






102 


1875 


9136 


- 34 






113 


2083 


10150 


- 64 






120 


2213 


10783 




-144 




130 


2396 


11675 


-144 






134 


2465 


12011 




-204 




142 


2604 


12688 


-204 






148 


2727 


13288 






-351 


165 


3036 


14793 




-354 




170 


3125 


15288 


-354 






t INCREASING 
i DECREASING 


PRESSURE, * MAXIMUM 
PRESSURE 


PRESSURE 
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Figure 19 Total System Setup 






864 OB Signal Generator which has very good spectral purity. 
C5575 could be driven to provide a maxiroum focal pressure of 
about one and a half atmospheres without using an amplifier. 
The voltage was measured with a HP 400E analog AC voltmeter 
or the oscilloscope. 

The excitation signal was originally generated by a 
HP 3 3 14 A Function Generator in the sweeping mode. This method 
produced bubble detection in most cases, but proved to be 
somewhat ambiguous in bubble size estimation. With the 
sweeping function generator, any given bubble's resonance 
frequency is broadcast once in the sweep interval. However, 
by broadcasting band limited white noise every frequency is 
present (though with randomly fluctuating amplitude) , 
increasing the probability of detection. A HP 8904A 
Multifunction Synthesizer was used to generate the band 
limited white noise. An Amplifier Research model 50A15 
amplifier with a maximum power of 50 watts and a frequency 
range of 20 kHz to 15 MHz was used to amplify the noise. A 
passive low pass RC filter was used to eliminate frequency 
components above 500 kHz. This signal was monitored 
periodically on the spectrum analyzer to verify proper 
frequency coverage. An oscilloscope was used to verify the 
overall noisy quality of the signal. The exciter was driven 
at high power to ensure enough energy was present at each 
frequency to excite any bubble in the size range of interest 
into large amplitude resonant oscillations. 
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2. Receiving Equipment 

The receiving system consists of a common branch 
supplying three independent detection devices. The receiving 
transducer converts acoustic signals to electrical signals. 
This signal passes through a DC blocking capacitor and an 
amplifier with a constant gain of 27.7 dB over the range of 
interest. Three paths diverge from this point. The first 
goes directly to an oscilloscope. The second into the 1 
megaohm input to the HP 3585A Spectrum Analyzer. The third 
path leads to a HP 3561 Dynamic Signal Analyzer via a 
frequency mixer. This third path was often disconnected due 
to certain limitations of the mixer. The functions of these 
detection devices are discussed in this section and the 
results chapter. 

The HP 3585A Spectrum Analyzer is the heart of the 
development version of this bubble detector. The oscilloscope 
was useful, but not capable of yielding quantitative 
information due to the short duration of the signal and the 
overwhelming magnitude of the scattered carrier compared to 
the sum frequency signal. A frequency mixer was employed to 
take the high frequency sum and difference signals down to 
base band. Then the HP 3 561 Dynamic Signal Analyzer could 
search the first 100 kHz to detect the larger, faster bubbles. 
This setup was both useful and frustrating. Neither of these 
instruments are continuously attentive to all frequencies. 
Bubbles were missed while these machines were busy with 
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calculations or searching other frequencies. Methods were 
hypothesized to correct this, but not implemented. These 
methods will be discussed with the conclusions. 

The oscilloscope would display the real time received 
signal. Without bubbles in the sample volume, this display 
was a sine wave at the carrier frequency. During a bubble 
event, this display would increase dramatically in amplitude 
as the bubble would scatter the exciter and carrier, as well 
as, the sum and difference frequencies toward the receiver. 
A sharp increase in the oscilloscope trace amplitude would be 
interpreted as a scattering event. It was not possible to 
tell what type of scatterer caused the event due to the speed 
of the fluctuation. This was used as an indication that a 
bubble actually passed through the sample volume. 

The HP 3585A uses a frequency sweeping receiver. The 
receiver is sensitive to a narrow bandwidth that sweeps across 
the frequency range specified. The fastest sweep time is 0.2 
seconds. A constant signal will be recorded each sweep. An 
intermittent signal is recorded only if the sweep is covering 
that frequency coincidentally. By adjusting the frequency 
range to sweep both the sum and difference frequency, each 
bubble has two chances to be detected each sweep. Hence the 
frequency range for bubble detection was 2.25 MHz to 3.05 MHz. 
This permitted detecting bubbles of up to 400 kHz resonance 
frequency. 
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The HP 3561 uses a fast Fourier transform (FFT) 
technique. It records an input signal of at least 4 
milliseconds, and then computes and displays the frequency 
spectrum of that signal via a FFT. It does not record data 
during the computation. An alternative to this method is to 
take 40 consecutive records of 4 milliseconds, then analyze 
and display the whole batch. This method is both better and 
worse. A greater period of time is covered during the 
consecutive records, but a greater period is left uncovered 
during the computations. 

The mixer required for the use of the HP 3561 has two 
distinct disadvantages. First, it distorts the received 
signal mildly in the region of the sum and difference 
frequencies as seen on the HP 3585A. Second, its own power 
limitations limit the power that can be simultaneously 
supplied to the carrier. Of these two, the distortion is the 
more serious problem. 

E. MODE OF OPERATION 

Many experiments were performed during the development of 
this system. Two major tests comprised the validation of the 
method to detect and size bubbles. The first involved timing 
a single bubble as it rose through a known distance. The 
resonance frequency corresponding to the rise time was then 
compared with the resonance frequency corresponding to sum and 
difference frequencies indicated on the HP 3585A. The purpose 
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of this test was to show that the sum or difference 



frequencies were being detected correctly, and that the sum 
or difference frequency gave an accurate and consistent 
indication of the bubble's resonance frequency. 

The other major validation test involved sending a single 
bubble through the sample volume and comparing the outputs of 
the HP 3585A, HP 3561 and oscilloscope. This test was 
designed to show that, at least for bubbles having resonance 
frequencies less than 100 kHz, the spectrum analyzer gave 
results consistent with those obtained from the HP 3561 Signal 
Analyzer. The two analyzers process signals differently. The 
spectrum analyzer samples a given frequency component only 
once per sweep. The signal analyzer digitizes a complete 
window of data and then performs a Fourier transform on that 
data. In that sense, the signal analyzer processes all 
frequency components simultaneously. 

The two validation tests had one significant limitation. 
Bubbles smaller than 30 urn radius cannot be easily seen by the 
naked eye to be timed and the HP 3561 cannot detect bubbles 
with resonance frequencies greater than 100 kHz. Therefore, 
resonance frequency determination of bubbles smaller than 30 
/xm radius (which also have resonance frequencies greater than 
100 kHz) cannot be verified by either method. In order to 
partially validate the system at the higher resonance 
frequencies, coincidence was used. If a very small radius 
bubble was detected by the system after the bubble generator 
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was momentarily energized, but no bubble was seen, then the 
system would be believed. This is not validation, but it is 
reasonable extrapolation. 

Other tests were performed with multiple bubbles in order 
to determine the systems response in this situation. Bubbles 
were readily detected, but the system display saturated 
quickly. This effect will be explained in the results 
section. 
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IV 



RESULTS 



The results of the tests discussed in the previous section 
are presented in this chapter. However, a brief discussion 
of the output display of the spectrum analyzer is presented 
first. The results of the single bubble detection tests, 
including data from the rise time resonance frequency 
verifications, follow. Tests of the dependence of the output 
signal level on the exciter and carrier pressure levels are 
then considered. Finally, the results of injecting solid 
particles into the sound field are discussed. 

A. TYPICAL OUTPUT DISPLAY 

The HP 3585A is the primary detection device for this 
system. This spectrum analyzer sweeps its frequency window 
through the range specified by the operator and displays a 
graph of amplitude versus frequency. Since the signal from 
a passing bubble is a fast acting transient phenomena, the 
most useful display mode for the purposes of this thesis is 
the "maximum hold" mode. In this mode the current display 
reflects the highest amplitude that has been present in a 
given frequency bin since the start of data acquisition. The 
use of this mode not only results in the display of the 
maximum signal detected but also the maximum noise. 
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In order to recognize a bubble detection, the display must 
be understood. Without bubbles, a small portion of the 
exciter and carrier signals are detected by the receiver. 
Figure 20 is a plot of the HP 3585A display of the received 
signal under this case. Figure 21 is the same display after 
the mixer has been added to the receiver circuitry. Notice 
that the mixer has both increased the noise in the frequency 
range 2.25 MHz to 3.05 MHz and reduced the 2.65 MHz carrier 
signal. Notice also the clean spectral nature of the carrier 
and the variations in the exciter level (0 4 00 kHz) . The 
high amplitude signal from 4 to 5 MHz in Figure 21 is higher 
order modulation of the carrier signal with the received 
signal. This noise is excusable because it is not in the 
bubble detection frequency range. 

In order to detect rapidly rising bubbles, the fastest 
sweep over the frequency region of interest was needed. This 
resulted in the 800 kHz span centered at 2.65 MHz being 
covered in a sweep time of 0.2 seconds. Any bubble with 
resonance frequency, in the range 3 0 to 4 00 kHz would have 
both a sum (2.65 MHz + f^) and difference (2.65 MHz - f^) in 
this frequency band. Figure 22 shows the display of this 
region as it appears without bubbles (and without the mixer) . 
Figure 23 shows the same frequency span after copious bubbles 
were allowed to rise through the sample volume for five 
minutes. A comparison of the two figures shows that the 
carrier level is at least 27 dB higher in Figure 23, 
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Figure 20 HP 3585A Display without Bubbles, Mixer 
Disconnected 
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Figure 21 HP 3585A Display with Mixer Connected, No Bubbles 
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Figure 22 HP 3585A Display of the Detection Band, without 
Bubbles, Mixer Disconnected (Carrier Focal Pressure 
= 5070 Pa, Exciter at Maximum Power 
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Figure 23 HP 3585A Detection Band Display with Bubbles, Mixer 
Disconnected (Carrier Focal Pressure = 5070 Pa, 
Exciter at Maximum Power) 




indicating that it was scattered significantly towards the 
receiver by the bubbles. Sum and difference frequencies were 
recorded for many bubbles of various sizes with a signal to 
noise ratio approaching 15 dB. Figure 24 is a superposition 
of the same type of plots for a higher carrier pressure. The 
bottom curve is without bubbles and without the mixer. The 
top curve is with bubbles but not the mixer. The middle curve 
has both bubbles and the mixer. As shown by the top curve, 
sum and difference frequencies are detected with a 10 to 20 
dB signal to noise ratio. These results are strong evidence 
which support the bubble detection goal. 

Figures 23 and 24 also illustrate system saturation. When 
the spectrum analyzer is in the maximum hold mode, the 
signature of a bubble, once detected, is maintained on the 
display. This signal is now noise to future detections and 
the system becomes less sensitive to bubbles of size similar 
to those it has already detected. This saturation effect is 
not a problem for these experiments since they involved only 
one bubble at a time. However, it indicates that this 
processing arrangement is not well suited for multiple 
detections . 

B. SINGLE BUBBLE DETECTION AND RESONANCE FREQUENCY 

DETERMINATION 

Figures 25, 26 and 27 show three variations of single 
bubble detection and resonance frequency determination. In 
each case, both the sum and difference frequencies were 
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Top: Yes Disconnected 26800 Pa 
Middle: Yes Connected 12800 Pa 
Bottom: No Disconnected 26800 Pa 
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Figure 25 Detection and Resonance Frequency Determination of 
a 55.3 kHz (by Rise Time) bubble. Peak Signal is 
at 56 Khz. Carrier Focal Pressure = 102400 Pa, 

Exciter at Maximum Power 
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Figure 26 Detection and Resonance Frequency Determination of 
a 62.4 KHz (by Rise Time) Bubble. Peak Signal is 
at 65.6 kHz. Carrier Focal Pressure = 102400 Pa, 
Exciter at Maximum Power 
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Figure 27 Detection and Resonance Frequency Determination of 
a Bubble too small to verify with the rise time 
method. Peak signal is at 100 kHz. Carrier Focal 
Pressure = 102400 Pa, Exciter at Maximum Power 



detected. Notice that the signal to noise level is a little 
different in each case. This difference is most likely due 
to the wide variety of carrier pressure available in the 
sample volume. 

The frequency of highest amplitude, other than the 
carrier, corresponds to the bubble's resonance frequency. 
Using the rise time method to independently determine 
resonance frequency verified the accuracy of the frequency 
determination. A large number of such comparisons were made. 
Figure 28 shows the agreement between the two methods. 
Perfect agreement would correspond to all of the data points 
falling on the line having a slope equal to one. The two 
methods agree to within a few kilohertz in most cases. A 
major source of error lies in timing the tiny rising bubbles. 

While performing these single bubble experiments it was 
noted that not every bubble was detected. There were some 
cases where bubbles passed through the sample volume 
sufficiently to show a transient on the oscilloscope but did 
not result in sum or difference frequency peaks on the HP 
3585A display. Since the oscilloscope verified the presence 
of a received signal, the HP 3585A was apparently inattentive 
to the frequencies of interest. The absence of these 
frequency peaks could have been due to: (1) the HP 3585A 
performing an internal calibration procedure that prohibited 
sampling while the bubble was in the sample volume, (2) the 
HP 3585A sweeping the wrong frequencies during a fast bubble 
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Rise Time Resonance Frequency (kHz) 




Figure 28 Resonance Frequency Comparison 
Rise Time versus HP 3585A 
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transient, (3) the bubble passing through a portion of the 
sample volume having too low a signal to noise level to permit 
detection, and (4) too much energy being scattered into the 
receiver causing the HP 3585A to overload. Some credence is 
given to all of these possiblities. Some displays would 
detect the sum or difference frequency but not both. Only 
alternative number (2) above makes sense in this case, since 
both signals should be present in approximately the same 
magnitude at the same time. 

The HP 3561 uses an FFT technique to analyze the spectral 
nature of the signal. However, the use of this instrument did 
not improve the probability of detection for several reasons. 
First, there is significant time over which the instrument is 
not sampling during computations, calibration and display. 
Second, the instrument's frequency range is limited to 0 to 
100 kHz. Hence, it could never detect a 100 to 400 kHz 
bubble. Finally, it requires use of the mixer which 
inherently lowers the signal to noise ratio. Some low power 
detections might be lost in the noise. 

However, the HP 3561 was useful. First, it verified that 
FFT techniques can be employed in this problem. Second, it 
provided a different method to verify the accuracy of the 
resonance frequency determination. Figure 29 shows the 
display of the HP 3561 for the case without bubbles. This 
shows the noise floor as a function of frequency. Since the 
mixer demodulates the signal from a 2.65 MHz carrier, the 
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Figure 29 HP 3561 Display without Bubbles 



frequency scale labeled 0 Hz to 100 kHz actually represents 
2.65 to 2.75 MHz. 

Figure 30 shows the detection and resonance frequency 
determination of a single bubble by the HP 3561. Using the 
time buffer mode, the HP 3561 continuously recorded 40 data 
records prior to processing. Ten of these data records were 
analyzed and plotted here. The first three and last six of 
the ten records show no hint of a bubble. Record four, on the 
other hand, strongly indicates the presence of a 40.25 kHz 
bubble. This figure indicates the transient nature of the 
bubble signals. The HP 3585A simultaneously detected the same 
bubble. Figure 31 is a copy of the HP 3585A display for this 
run. It shows the difference frequency strongly at 2.61 MHz, 
which represents a resonance frequency of 40.00 kHz. The sum 
frecjuency did not show up due to the high speed of this buble 
transient. In other words, the signal was gone before the 20 
milliseconds needed for the sweep to advance to the sum 
frequency had elapsed. The HP 3561 bears this out since each 
record represents only 4 milliseconds of data. Although two 
peaks of approximately the same amplitude are shown in Figure 
30 the one at 40.25 kHz is actually much larger. It appears 
to be the same amplitude because, it exceeded the maximum 
scale and was truncated. The level -87.13 dBV is, however, 
accurate . 
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Figure 30 HP 3561 Display of 10 Data Records’ Frequency 
Content. The Fourth Trace Shows a 40.25 kHz Bubble 
Detection . 
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Figure 31 HP 3585A Display Showing Detection of the Same 
bubble as Figure 30. Peak Signal is at 40.0 kHz 



PRESSURE 



C. VERIFICATION OF DUAL FREQUENCY METHOD 

DEPENDENCIES 

In order to confirm that the detected signals are in fact 
produced by the mechanism described by the dual frequency 
method, the peak pressure amplitudes of the detected signals 
were compared to the predictions of equation 2. The 
calculations were performed for many different bubble sizes 
using the transducer calibrations to determine the maximum 
carrier pressure, the exciter pressure at the bubble resonance 
frequency and the pressure amplitude of the sum or difference 
frequency. Figure 32 is a comparison of the actual received 
level and maximum calculated sum or difference pressure at the 
receiver face divided by the product of carrier and exciter 
pressure. The expected values were given earlier in Table 4. 
None of the actual values exceeded the maximum expected 
values. The lower actual values are attributable to the fact 
that most bubbles encountered less than maximum carrier 
pressure. 

According to the theory of dual frequency bubble 
detection, the amplitude of the return signal depends on the 
pressure amplitudes of the carrier and excitation at the 
location of the bubble. In order to get an idea of what drive 
pressures were required, the exciter pressure was varied over 
a range of 7.36 dB with a constant carrier pressure amplitude 
of 102 kPa. Since the carrier pressure varies dramatically 
over short distances and the excitation pressure varies wildly 
with frequency, it would take an enormous amount of data to 
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Figure 32 Comparison of Actual Pressure Relationship 
with Maximum Expected 
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verify equation 2 with this implementation. In spite of this, 
Figure 33 shows that the signal level varies with the exciter 
pressure. The exciter pressure is lowest for frequencies of 
about 140 kHz. Therefore, 140 kHz bubbles are detected with 
the weakest signal. The main reason for the curves not being 
flat is the variations of excitation pressure with frequency. 
Another factor is that smaller bubbles naturally have smaller 
signals for the same excitation. 

Figure 34 shows a similar family of curves. This time the 
exciter was held constant at the maximum value. The carrier 
pressure was varied by a factor of 2 (6dB) . Again, the bubble 
signal decreased as carrier pressure decreased. It is clear 
that for better signal to noise ratios, the carrier and 
exciter pressures should be as high as possible. 

D. DISCRIMINATION OF BUBBLES FROM SOLID PARTICLES 

Dry grains of sand, ranging in size from 1 ^m to 1 mm 
diameter were dropped into the water over the sample region. 
The grains rarely entered the sample volume, but when one did 
it almost always produced a signal similar to that produced 
by bubbles. These signals were attributed to the presence of 
entrained gas pockets, some of which were visible, on the dry 
sand. It was reasoned that wet sand would not produce the 
same type of signals. Therefore, pre-wetted sand was more 
carefully injected into the sample volume. Only occasionally 
would an extraneous bubble indication be noted under these 
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Received Voltage Level (dBv) 




Figure 33 Received Signal Variations with Changing 
Exciter Pressure 
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Received Voltage Level (dBv) 




Figure 34 Received Signal Variation with Changing 
Carrier Pressure 
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circumstances. This was thought to be due to an infrequent 
gas bubble in the slurry. The most common response of the 
system to a solid is an increase in the received carrier and 
exciter signals. As long as the pressure of these signals is 
less than that required for medium non-linearities, such as 
the parametric effect, no sum or difference frequencies should 
be detected. 

Hollow glass beads were also injected into the sample 
volume. The response was identical to that for gas bubbles. 
This response semed reasonable since the glass beads were at 
least 90 percent gas by volume. They floated rapidly to the 
surface. The results of these measurements indicate that the 
dual frequency method is capable of detecting bubbles as well 
as other partially gaseous matter. However, it is insensitive 
to non-gaseous scatterers. 

E. DISCUSSION 

Up to this point, many facts and figures have been 
presented, but the full implications have not been discussed. 
This section will review the evidence presented for this 
implementation of the dual frequency method for detecting and 
sizing bubbles. 

The first requirement was to detect bubbles. This goal 
clearly was accomplished. Bubbles were detected through the 
range of interest, though the majority of the events were less 
than 100 kHz. Data to support this was taken with single 
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bubbles and lots of bubbles, free rising bubbles and convected 
bubbles, varied excitation and carrier pressure levels. By 
doing this, not only were bubbles detected, but some 
limitations of bubble detectability were discovered. These 
will be highlighted later in this section. 

The ability to directly determine the resonance frequency 
of a bubble has been shown. Two independent signal processors 
were used to show that this technique was not based on a 
peculiarity of one type of processor. The resonance frequency 
determination follows the simple rule that the bubble has its 
maximum response at resonance. Therefore, the maximum 
received signal in the detection band corresponds to the 
resonance frequency. The rise time method of resonance 
frequency estimation supported the dual frequency method 
results. The agreement in the results of the two processors 
and the rise time method was best when detected with high 
signal to noise. This was not surprising since low signal to 
noise was often due to poor bubble signal timing with respect 
to the frequency sweep for the HP 3585A or the signal data 
record for the HP 3561. Strong signal to noise indicated that 
the bubble was properly detected while in the processor's 
analysis window. In these cases, this dual frequency method 
determines the resonance frequency to within the signal 
processor's analysis error. 

The tests with solid particles indicated that this method 
can determine the resonance frequency of solid borne or 
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roembrane enclosed gas volumes. The accuracy of these 
measurements has not been shown, but the presence of a signal 
very similar to that of a pure bubble for these gas volumes 
with solids suggests the same ability to correctly determine 
resonance frequencies. A solid particle without gas was shown 
to lack the sum and difference frequency signals. Therefore 
non-gaseous scatterers are easily distinguishable from various 
types of gaseous scatterers having resonance frequencies in 
the range of interest. This natural ability to classify a 
bubble by its acoustic resonance parameters has far reaching 
implications for the study of bubbles in both clean and dirty 
environments . 

During the bubble experiments, a wide variety of signal 
to noise levels were noted. Some bubbles transited the sample 
volume without any indication of its presence on the 
processor's display. Others were detected with more than 20 
dB of signal to noise. Much of this was probabilistic since 
the bubble's path through the sample volume and signal timing 
to the signal processor were not controllable. This random 
positioning and timing meant that any single run with any 
single bubble could not be directly related to other runs. 
To compare one test setup with another many data runs were 
taken. The maximum signal to noise cases were considered to 
have had the most optimal sample volume positioning and 
timing. Thus, the differences between these runs would be due 
to the actual parameter changes between test setups. An 
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example of a test setup change is the exciter pressure 
variations previously discussed and shown in Figure 33. To 
verify the relationship of bubble signal pressure to exciter 
pressure many bubble runs were made at each pressure. The 
maximum bubble signal at each pressure was considered to be 
optimum in position and timing. 

The sum frequency pressure relationship given by equation 
2 was tested by doing these comparisons. To establish a 
threshold of detection the whole system must be understood. 
For this system, the detection band noise level was about -100 
dBv. To obtain 10 dB signal to noise, a signal of -90 dBv is 
needed at the instrument. Subtracting the 27.7 dB amplifier 
gain yields a necessary -117.7 dBv out of the receiver. It 
is easy to calculate the pressure required at the receiver 
face to generate this voltage. Using equation 2, the required 
excitation pressure can be found for each frequency given the 
carrier pressure. Table 14 shows the results of these 
calculations . 

The pressure produced by the exciter was well above the 
calculated values of Table 14 . The carrier pressure was also 
higher than the value assumed for this table. These higher 
pressures produced the data with 20 dB signal to noise. The 
weaker signals were due to the losses induced by poor 
positioning and timing. The numbers in Table 14 are valid for 
any system since they are simply pressure relationships. The 
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TABLE 14 


EXCITER PRESSURE REQUIRED 
GIVEN r = 7.7 cm, CARRIER 


FOR 10 dB SIGNAL TO 
PRESSURE = 3 X 10^ 


NOISE 

Pa 


FREQ. 

fMHz) 


P+ at Receiver for 
-117.7 dBv fPafoeak^l 


Exciter 
Pressure 
f Pa ^oeak ) ) 


FREQ. 

^kHz) 


2.70 


0.0197 




29.6 


50 


2.75 


0.0194 




68.0 


100 


2.80 


0.0195 




115 


150 


2.85 


0.0191 




184 


200 


2.90 


0.0191 




254 


250 


2.95 


0.0190 




274 


300 


3.00 


0.0189 




319 


350 


3.05 


0.0190 




374 


400 



113 



receiver output voltage would, of course, vary with different 
receivers. This makes this analysis specific to this system. 

If a system could be designed to produce these exciter and 
carrier pressures and receiver coverage over a broad sample 
volume, it would detect all bubbles with about 10 dB signal 
to noise ratio. This assumes that a larger sample volume 
would alleviate the critical timing problem since the bubble 
would remain in the sample volume longer. 
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V. CONCLUSIONS 



The goals of this bubble detection system have been 
accomplished. The system has demonstrated the ability to 
detect bubbles with resonance frequencies from 30 to 400 kHz 
when they entered the sample volume. The resonance 
frequencies of these bubbles were accurately determined by the 
system. The ability to distinguish between bubbles and nonr 
gaseous particles was shown. Parameters for creating a system 
with a high probability of bubble detection were accumulated. 
The ability to determine the resonance frequency of scatterers 
containing a small gas volume was also noted. 

With the achievement of these goals, the dual frequency 
method for bubble detection and resonance frequency 
determination becomes a powerful tool. It can be used to 
analyze bubble populations or to study the dynamics of single 
bubbles under the influence of stress. Studies may be done 
using this technique to study the effect of varying a fluid's 
viscosity or surface tension on a bubbles damping constant. 
There are a number of areas for application of this technique. 
The next section recommends items to consider to improve the 
performance of a dual frequency method bubble detector. 
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A. RECOMMENDED IMPROVEMENTS FOR A CONTINUOUS BUBBLE DETECTION 
SYSTEM 

Many items have been mentioned that effect the chances of 



detecting a bubble with this 
be made to achieve an ideal 
areas of concern. 

* Sample Volume 

* Bubble 

* Receiving Equipment 



* Environment 



system. Some improvements must 
system. Here is a summary of 



Concerns 

- too small 

- too variable 

- too weak 

- path and time in sample volume 

- needs to sample continuously 

- needs large dynamic range (to 
detect small sum and 
difference frequency signals 
in the presence of large 
scattered carrier and exciter 
signals) 

- data output management 

- noise level in detection band 

- number of bubbles: saturation 

effect and signal extinction 

- vibrations: disturb sample 

volume alignment 

- flow rate: speed of bubble 

through sample volume 

- quantity of suspended solid 
scatterers 



1. Sample Volvune Related Improvements 

A bubble must pass through the sample volume to be 
detected. A large sample volume permits easier coverage of 
reasonable volumes. The smaller the sample volume is with 
respect to the total system, the lower the probability of 
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detection. The implementation tested here had virtually the 
smallest possible sample volume. 

Variations in the carrier and excitation pressures and 
receiver efficiency within the sample volume allow for highly 
variable received signals when a bubble is present. If the 
lowest possible received signal is less than the noise, 
detection opportunities will be missed. Driving the exciter 
with band limited noise in the frequency range of interest 
provides broad frequency coverage, with only minor temporal 
variations. The drive current and exciter source strength are 
functions of frequency. The drive current should be 
manipulated to excite bubbles equitably, to give each size 
bubble the energy to obtain the same signal to noise. 

The receiver's efficiency may vary with frequency and 
bubble location. These variations are strictly functions of 
the transducer. The transducer should be designed or selected 
to minimize these variations in the frequency range and 
spatial region of interest. The carrier transducer is 
monofrequency. Hence, only stability of output power and beam 
pattern are of concern. The overall goal is to obtain a good 
sized volume of constant and adequate, but not excessive, 
acoustic pressure. If this is accomplished, many of the 
problems associated with a small sample volume disappear. 
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2. Bubble Related Concerns 



The path of travel varies from bubble to bubble, due 
to a variety of forces acting on the bubble. Gravity, 
buoyancy, drag and radiation forces are the key players. Of 
these, radiation force may be the worst. This acoustic 
phenomena pushes the bubble when the acoustic pressure is 
high. The resulting motion drives the bubble toward regions 
of lower carrier and excitation pressure. Keeping the 
pressure just below the level that causes significant motion 
is important to solving this problem. Having a relatively 
large sample volume would help regardless of the cause of the 
bubbles motion. 

The amount of time a bubble spends in the sample 
volume impacts its detectability. It is a function of the 
sample volume, bubble speed and path. These items need to be 
manipulated to achieve an ideal length of time in the sample 
volume. Too little time results in missed detections. Too 
much time will limit the rate at which the system can 
distinguish between individual bubbles. 

3 . Receiving Equipment Improvements 

The receiver must be alert to detect an event. The 
equipment assembled for this experiment was useful, but not 
ideal for the job. The received signal needs to be 
continuously monitored over the frequency range of interest. 
Only the sum or difference band is actually needed, but both 
may be wanted for verification. The data needs to be 
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presented in a format in which new detections are not effected 
by previous ones. A continuous hard copy spectrogram is a 
good way to handle large amounts of data. 

During a bubble detection, many frequency signals are 
received. Unfortunately, those of the most interest are of 
low level. Large dynamic range is needed in the receiver. 
This allows it to measure the low level signals without 
overloading due to the high level scattered signals. 
Filtering the return signal can help this problem. Band pass 
filtering to accept only the detection band followed by notch 
filtering to eliminate the carrier would eliminate the 
unwanted signals' overbearance. 

Correcting the overall alertness and relative dynamic 
range of the receiving equipment would be a major improvement 
in this system. This would require a system designed for this 
purpose. Most off the shelf equipment does not have the 
frequency range or data volume capability needed for a true 
bubble detection system. 

4. Environment Related Concerns 

The environment of the test apparatus has a large 
impact on its effectiveness. Vibrations change the bubbles 
path, as well as, the sample volume alignment. The quantity 
of scatterers throttles the amount of power available to 
excite a bubble and effects the propagation of the bubble's 
signal to the receiver. The ambient noise in the detection 
band is a factor in detectability. The environment was not 



119 



a major problem in this laboratory experiment. It has the 
possibility of causing problems, but they are minor compared 
to the previous three categories. 
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APPENDIX 



EASY REFERENCE TABLE 



BUBBLE RES . FREQ 

RADIUS, a fum) fkHz> 



DAMPING 1 

CONSTANT Pi P? Pa 



TIME TO 
RISE 1"(SEC^ 



100 


31.4 


0.075 


70 


44.6 


0.085 


50 


62 . 1 


0.095 


40 


77.3 


0.100 


30 


102.6 


0.115 


20 


152.9 


0.13 


10 


303.8 


0.16 


9 


337.7 


0.16 


8 


380.4 


0.16 


7 


435.8 


0.16 


6 


510.8 


0.16 


5 


618.0 


0.16 



4.3 


X 


10 '® 


1.44 


2.7 


X 


10 '® 


2.55 


1.74 


X 


10 '® 


4 . 61 


1.34 


X 


10 ‘® 


6.95 


8.79 


X 


10 ’’ 


12.0 


5.25 


X 


10 '’ 


26.50 


2.16 


X 


10 ’’ 


105 


1.94 


X 


10 "’ 


129 


1.72 


X 


10 ’’ 


163 


1.50 


X 


10 ’’ 


213 


1.28 


X 


10 ’’ 


— 


1.05 


X 


10 ’’ 


416 



EXCITER 



FREQ 


PpEAK 


fkHz) 


(Pa) 


50 


115 


100 


261 


150 


244 


200 


370 


250 


515 


3O0 


980 


350 


1000 


400 


1190 



CARRIER 



Vpp 


PpEAK 


( Pa ) 


0.052 


707 


0.165 


2242 


0.41 


5570.0 


0.55 


7472.0 


0.94 


12800 . 


1.10 


14940 . 


1.40 


18250 . 


1.98 


26820 . 


2.08 


28260 . 


2.55 


34640 . 


3 . 90 


51800 . 


4.0 


47080 . 


4.05 


55020 . 


5.0 


66580 . 


5.2 


70650 . 


6.0 


81540 . 


7.0 


88080 . 


7.67 


102400 . 


8.0 


105300 . 


8.2 


111400 . 


9.0 


115300 . 


10.0 


133200 . 


10.35 


140600 . 


11.0 


148900 . 


12.9 


175300 . 



RECEIVER 

FREQ M (V/Pa) 



fMHz) 


2 . 25 


0.0000982 


2.3 


0.0000969 


2 .35 


0.0000955 


2.4 


0.0000951 


2 .45 


0.0000939 


2.5 


0.0000936 


2 . 55 


0.0000932 


2.6 


0.0000934 


2 . 65 


0.0000932 


2 . 7 


0.0000937 


2.75 


0.0000948 


2 . 8 


0.0000947 


2.85 


0.0000967 


2.9 


0.0000966 


2.95 


0.0000982 


3 . 0 


0.0000974 


3.05 


0.0000972 


3.1 


0.0000966 
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